First-principles Calculation of the Single Impurity Surface Kondo Resonance 
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We have performed first-principles calculation of the surface and bulk wavefunctions of the Cu(lll) 
surface and their hybridization energies to a Co adatom, including the potential scattering from the 
Co. By analyzing the calculated hybridization energies, we found the bulk states dominate the 
contribution to the Kondo temperature, in agreement with recent experiments. Furthermore, we 
also calculate the tunneling conductance of a scanning tunneling microscope (STM) and compare 
our results with recent experiments of Co impurities in the Cu(lll) surface. Good quantitative 
agreement is found at short parallel impurity-tip distances ( < 6 A). Our results indicate the need 
for a new formulation of the problem at larger distances. 

PACS numbers: 68.37.Ef, 72.15.Qm, 72.10.Fk 
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Scanning tunneling microscopy (STM) has become one 
of the most basic tools for the manipulation of mat- 
ter at the atomic scale. Although this experimental 
technique has reached maturity, the detailed theoretical 
understanding of experimental data is still incomplete 
and/or contradictory. One of the most famous exam- 
ples of atomic manipulation is associated with the surface 
Kondo effect observed when transition metal ions (such 
as Co) are placed on a metallic surface (such as Cu(lll)) 
P| . The surface Kondo effect is the basis for the observa- 
tion of surprising phenomena such as quantum mirages 
0, and has attracted a lot of attention and interest in 
the last few years. The current understanding of these 
observations is based on the assumption that only sur- 
face states of Cu(lll) are involved in the scattering of 
electron waves by the Co adatoms [10113 Neverthe- 
less, recent experiments with Co atoms on the Cu(fOO) 
surface (that does not have any surface state) [6j, or in 
Cu(lll) but close to atomic surface steps (that affect the 
surface states) |j| have indicated that bulk (not surface) 
states are behind the surface Kondo effect. Meanwhile, 
in contrast, the growing theoretical literature in the sub- 
ject is heavily concentrated on the surface states alone. 
In this paper, we use first-principles methods that clearly 
show that the bulk states are behind the surface Kondo 
effect, in agreement with these experiments. In the light 
of these results, theoretical approaches based on surface 
states alone have to be reconsidered. 

The surface Kondo effect without the STM probe is 
described by the Anderson impurity model: 

H s = ^2 £kcl a c ka + ^2(t ka cl a c aa + H.c.) 

k<7 ktr 

+ ^2 e aci c7 c aa + Un a - l n al , (I) 

a 

where ek is the energy of the substrate electrons, c fc 
(cku) creation (annihilation) for substrate electrons with 
momentum k and spin a, e a is the energy of the adatom 
d-orbital, (c aa ) is the creation (annihilation) of 



adatom electrons, and U is the Coulomb energy for dou- 
ble occupancy of the d-orbital (n acr is the number oper- 
ator). The hybridization energy between the substrate 
and adatom is written in terms of matrix elements be- 
tween their wavefunctions. 



t ka = f d 3 r^(r)(H + V a (r))^ a (r), 



(2) 



where Hq refers to the bare metal surface, V a (r) is the 
adatom scattering potential, V'k(r) is the substrate wave- 
function scattered by the adatom potential, and ip a (r) 
the adatom 3d orbital. Most work on the surface Kondo 
effect ja, |2( follows Anderson's original idea [iJJ] in deal- 
ing with the hybridization matrix element t^ a , namely, 
treating it as a phenomenological parameter to fit exper- 
iments. However, in trying to understand the STM ex- 
periments, and especially the role played by the surface 
and bulk states, these matrix elements cannot be simply 
taken as phenomenological parameters since one would 
not be able to separate the contributions coming from 
the bulk and the surface states of the substrate. Thus, it 
is essential to perform microscopic calculations of these 
matrix elements starting from the electronic wavefunc- 
tions. 

Microscopic calculations of these matrix elements using 
the nearly free-electron model (NFE) for Cu have been 
attempted recently [ill H^ |. The NFE approximation 
has its advantages in analyzing the momentum depen- 
dence of the hybridization energies and obtaining analyt- 
ical substrate wavefunctions. However, in using the NFE 
one needs to prescribe how the substrate crystal poten- 
tial joins its vacuum image counterpart, and the unique- 
ness of this potential prescription is questionable. Such 
a freedom in modeling potential makes the NFE method 
unreliable in obtaining the precise surface and bulk wave- 
functions. Moreover, we further notice that the substrate 
states in Q are those already scattered by the adatom 
potential rather than the bare-substrate states. The rela- 
tion between these adatom-scattered states and the bare 
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substrate states is given by: 

Et (°)t (o) 

+ EK'«l+^-), (3) 

kk'o- 

c/kk' = | dVv^^WKWv^W, (4) 

where the superscript (0) refers to the bare electronic 
states. The NFE does not consider the scattering poten- 
tial from the adatom and uses the bare-substrate states 
in the Anderson impurity model 

In the presence of the tip new terms have to be 
added to the Hamiltonian that describe the tunneling 
processes of tip-to-adatom and tip-to-substrate 0- The 
tip-to-adatom tunneling process is described by: H at = 
J2a- tap ( c L c p<t + H.c.) , where c pa (c pa ) annihilates (cre- 
ates) electrons at the tip and t ap is the hybridization 
energy between tip and adatom. The tip-to-substrate hy- 
bridization is given by: H st = J^ka- ^ P (cLr'V + H - c -) > 
where i kp is the hybridization energy between substrate 
(bulk or surface state) and adatom. The tip Hamilto- 
nian is simply: H t = e p c P a c pcn where e p is the energy 
of the tip electrons. The total Hamiltonian of the tip- 
substrate-adsorbate system is: H = H s +H t + H at + H st . 
For the hybridization energies that involve the STM 
tip, we follow Plihal and Gadzuk [9| and approximate 
t pa oc ?/>*(R t ) with a = k, a. 

To correctly obtain the contributions of the surface 
Kondo resonance from the surface and bulk states, we 
perform first-principles calculations of surface and bulk 
wavefunctions on the Cu(lll) surface in the presence of 
the scattering potential, as well as their hybridization 
energies to the Co adatom. As the first step, we calcu- 
late the wavefunctions of a bare Cu(lll) surface. Such 
a surface is simulated by a super-cell of 21-layer slabs 
separated by 8 vacuum layers. We employ, in the frame- 
work of density functional theory, a self-consistent full- 
potential linearized augmented plane wave (FLAPW) 
method |l3|. with the exchange-correlation potential in 
the generalized gradient approximation (GGA) 14]. The 
interatomic distances within a Cu slab are determined by 
the bulk lattice constant a = 3.62 A, and the surface 
relaxations of Cu(lll) are neglected because of its close- 
packed structure. The calculated surface-state dispersion 
agrees very well with the experiments (Fig. left). 

In order to obtain the potential of a Co adatom, we per- 
form another FLAPW calculation, in the local spin den- 
sity approximation (LSDA) , of 7-layer Cu slabs separated 
by 8 vacuum layers, plus Co impurities 1-layer spacing 
high on top of the Cu surface layers (Fig. fright). The 
potential V a {r) appearing in both J2J) and (@J is taken 
to be the potential difference between the surface with 
Co and the clean crystal. Using the bare Cu(lll) states 
within the energy range \e — €f\ < 1 eV as basis, we per- 
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Figure 1: The left figure shows the surface-state dispersion 
from our calculation (curve) and the experiments [lf| (data 
points). The right figure shows the top view of the unit cell 
of Co on the Cu(lll) surface in our LSDA calculation of the 
present work. The Co atom is located at the center (grey 
circle). The larger and smaller circles are the first (surface) 
and second Cu layers, respectively. The individual Co atoms 
are separated by 8 A. 

form an exact matrix diagonalization on the right-hand 
side of and obtain the adatom scattered states. It is 
the 3d orbital of the Co adatom that actually participates 
in the Kondo resonance. The 3d orbital c aa appearing 
in is renormalized by the Cu-substrate potential but 
does not include hybridization with the substrate states. 
Our current first-principles approach cannot generate 3d 
orbitals satisfying both conditions. Instead we can per- 
form LSDA calculations for the 3d orbital of either a sin- 
gle Co atom or a Co atom on Cu(lll) with hybridization 
i a k included. The latter is not a good candidate to be 
used in Q because it adsorbs i k into itself and should 
in principle give i ak ~ 0. Thus we calculate the elec- 
tronic structure of a single Co atom using a relativistic 
atomic code and use its 3d 32 2_ r 2 orbital as the ip a (r) 
in The particular choice of the 3d orbital 3d^2_ r 2 
other than 3d xy or 3d x 2 is supported by Ref. |l7j . 

The theory of the surface Kondo resonance adopted 
in the present work closely follows Ref. 0, The 
broadening of the Co d-level is calculated directly from 
ika in © (without any adjustable parameters) 

A = ^|< kQ | 2 <5(e f -e k ), (5) 

k 

where £f is the Fermi energy, to be A = 0.18 eV. In fact, 
the STM-measured Kondo temperature for Co/Cu(lll) 
is Tk ~ 50X, which, from Tk — Dexp (— 7r|e a |/A), gives 
A ~ 0.2 eV if one uses well established values of the 
d level e a ~ 0.9 eV Q and the Cu band cutoff D ~ 
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5.5 eV [l9|. The contributions to A from the surface 
and bulk states are also investigated by considering A = 
A sur f + A bulk . We calculate from (JSJ the ratio of the 
surface-state contribution to the total <i-level broadening 
A sur f/A ~ 0.006. The ratio shows that the bulk states 
dominate the adatom-substrate hybridization energy of 
Co/Cu(lll). This result can be understood by the fact 
that the Co atom is still in one layer above the surface 
Cu layer, a crystal-like regime rather than the tunneling 
regime (that is, the Co adatom is strongly hybridized 
with the substrate). 

The STM differential conductance can be written as: 



dl 
dV 



dl 

dV 



a{R) 



\q(R)\ 2 -I + 2£Re [q(R)} 



1 



(6) 



where R is the parallel impurity-tip distance, £ = (eV + 
£a) I '{^bTk) is the dimensionless bias, e is a bias off-set 
due to the d-state energy, 



a(R) 



7T ^ ipkika<5(£F - Ck" 



is the amplitude of the resonance, and 

_ tpa + X)k ^pk^kaRe Gk 
Sk *pk*koIm Gk 



(7) 



(8) 



is the so-called Fano parameter where Gk = (cf — £k + 
irj)^ 1 is the bare substrate Green's function. In 10, 
the differential conductance with a subscript "0" refers 
to the background signal (proportional to the local den- 
sity of states of the substrate). Typically q(R) has been 
taken to be real in the experimental fits. However, 
when performing the first-principles calculation of the 
Cu(lll) surface, we found that q can carry an imag- 
inary part. The Bloch states of the conduction elec- 
trons can generally be written as V , nk(r) = Mnk(r)e~ 4k|1 r u 
with M„k(r) = \u n u(r) \ e^^ r \ NFE studies of the surface 
Kondo resonance js|, [ill treat u ra k m approximation 
such that $(r) is spatially independent, i.e., $(r) = $o 
is an overall phase. However, our first-principles calcu- 
lation shows that the function u„k can carry a spatially 
varying phase. Including $(r) in the calculation of the 
hybridization energies consequently gives the complex q. 

The STM tip is positioned about 5 A< z < 10 A above 
the surface in the usual spectroscopic tunneling condi- 
tions. However, when performing first-principles calcu- 
lations of the bare-Cu(lll) wavefunctions, we find that 
the wavefunctions undergo an oscillatory behavior rather 
than a smooth exponential decay beyond 4 A from the 
surface. It is known that this problem comes from GGA 
in the low density region because of large-scaled gradi- 
ents H3- The large gradients cause fluctuations in the 
exchange-correlation potential, which leads to fluctua- 
tions in wavefunctions. To resolve this problem we ex- 
trapolate the Cu-substrate wavefunctions by fitting their 
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Figure 2: Normalized Fano formula prefactors a(R)/a(0) at 
three tip heights Zt = 3.5 A, 10 A, 16 A (from the bottom to 
the top). The experimental data is also shown 0. 



values in 2 A< z < 3 A using the NFE wavefunctions (llj| . 
The NFE wavefunction is then used in the z > 5 A region 
to calculate the STM differential conductance. 

We analyze surface and bulk-state contributions to 
a(R), and find that the bulk states dominate. In Fig. [21 
we plot a(R) of three tip heights (Z t = 3.5 A, 10 A, 
16 A) as well as the experimental data The pref- 
actor a(R) at Z t — 3.5 A is calculated directly from the 
GGA wavefunctions rather than the extrapolated NFE 
wavefunctions. It is clear that a(R) moves towards the 
experimental data as the tip moves farther away from 
the surface. One can see that our first-principles calcu- 
lation at Z t — 16 A (without fitting parameters) agrees 
well with the experiment within R < 5 A but starts de- 
viating from the experiments as R increases and has a 
node around R = 7 A. Since the current first-principles 
approach in calculating the Cu(lll) electronic structure, 
GGA, is widely regarded as a highly accurate compu- 
tation scheme except for very low electron density, we 
suggest that further theoretical and experimental work is 
required to check this issue. In the bulk Kondo problem 
the conduction-impurity hybridization can be regarded 
as momentum- independent, and the Kondo Hamiltonian 
has an exact solution by Bethe ansatz 21] . However, the 
surface Kondo effect has a k-dependent ik a , and there 
is so far no field-theoretical approach that can treat it 
exactly. A possible solution is to use a computational 
scheme to compute the Anderson Hamiltonian of a sur- 
face Kondo effect. In the experimental aspect it should 
be pointed out that at large values of R the STM data is 
very noisy, and it is possible that the fitting is not unique. 

The calculation of a(R) requires only the Cu-substrate 
states at the Fermi energy (see Eq. J7J)) while the line- 
shape parameter q in JSJ) depends on the entire Cu band. 
We use the calculated Co-scattered Cu states within the 
energy range |e — ep\ < 0.9 eV to obtain the q(R) in 
Eq. (JHJ. Since the Fano lineshape parameter q defined 
in (JSJ) is in general complex, the experimentally fitted q 
based on the Fano formula of JHJ with lm[q] = can only 
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Figure 3: The absolute value (solid curve) and real part 
(dashed curve) of calculated complex q at Zt = 16 A as a 
function of R as compared with the experimental fit assum- 
ing Im[g] = 0. Inset: comparison of theoretical (light) and 
experimental (dark) STM dl /dV vs. V lineshapes at 7? = 0. 

be compared with our calculation qualitatively. We plot 
our calculated \q\, Re q, and the experimentally fitted q 
in Fig. |21 The inset is a direct comparison of between 
our calculated and the experimental STM lineshapes at 
R = 0, showing good agreement. One can see likewise 
that the q vs. R plot shows good agreement between our 
calculation and the experiments for R < 6 A and the dis- 
crepancy for R > 6 A is a consequence of the node of 
a(R) around R = 7 A. 

The effect of the potential scattering from the Co atom 
is also studied by calculating the rf-level broadening A 
and the STM lineshape prefactor a(R) without potential 
scattering, i.e., C/kk' = in ©. When potential scat- 
tering is neglected, the d-level broadening A reduces by 
8%, and its surface-state contribution slightly increases 
but is still small (A sur f/A ~ 0.025). We also found that 
without potential scattering the surface states dominate 
the contribution to the STM lineshape prefactor a(R). 
This is a drastic change from the potential scattering 
case where bulk states dominate. This change can be 
understood as follows: bulk states dominate the local 
density of states (LDOS) of the conduction electrons at 
the adatom site with and without potential scattering, 
which accounts for the d-level broadening. As LDOS of 
the bulk states decays away from the surface faster than 
the surface states, in the case of no potential scatter- 
ing the product of t a k (bulk dominated) and t p k (surface 
dominated) in Q turns out to be dominated by the sur- 
face states. When potential scattering is included, the 
bare-substrate bulk and surface states strongly mix with 
each other. As a result, LDOS of the bulk and surface 
states in the presence of potential scattering decay from 
the surface approximately in the same rate, and t a kt p k 
becomes dominated by the bulk states. 

In summary, we have calculated the hybridization ener- 
gies from the LDA wavefunctions of the Cu (111) surface 



and Co atom. The potential scattering of Cu conduction 
states from the Co adatom is included in determining the 
substrate-adatom hybridization energy. Our calculated 
d-level broadening from the above hybridization energy 
is in excellent agreement with the value determined from 
the STM-measured Kondo temperature. Our analysis of 
the contribution of the substrate-adatom hybridization 
energy from surface and bulk states shows that the bulk 
states dominate the Kondo temperature. We also calcu- 
lated the tunneling conductance of an STM tip for the 
Cu(lll) surface in the presence of a Co adatom. Our cal- 
culated conductance has quantitative agreement with the 
experiments at short parallel tip-adatom distance with- 
out any adjustable parameters. However, discrepancy 
appears as the parallel distance increases indicating that 
a new approach is required for this problem. 
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